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New data for the total cross section o"(e+e~ hadrons) in the charm and bottom threshold 
region are combined with an improved theoretical analysis, which includes recent four-loop 
calculations, to determine the short distance MS charm and bottom quark masses. The final 
result for the MS-masses, mc(3 GeV) = 0.986(13) GeV and mi,(10 GeV) = 3.609(25) GeV is 
consistent with but significantly more precise than a similar previous study. 



The strong coupling constant and the quark masses are the fundamental input parameters of 
the theory of strong interaction. Quark masses are an essential input for the evaluation of weak 
decay rates of heavy mesons and for quarkonium spectroscopy. Decay rates and branching ratios 
of a light Higgs boson, suggested by electroweak precision measurements, depend critically on the 
masses of the charm and bottom quarks, nic and m^. Last not least, confronting the predictions 
for these masses with experiment is an important task for all variants of Grand Unified Theories. 
To deduce the values in a consistent way from different experimental investigations and with 
utmost precision is thus a must for current phenomenology. 

A detailed analysis of rric and mh based on the ITEP sum rules^^ has been performed several 
years ag(P and lead to mdmc) = 1.304(27) GeV and mb{mb) = 4.191(51) GeV. During the 
past years new and more precise data for a{e~^e~ hadrons) have become available in the low 
energy region, in particular for the parameters of the charmonium and bottomonium resonances. 
Furthermore, the error in the strong coupling constant as{Mz) = 0.1189 it 0.0020 has been 
reduced. Last not least, the vacuum pol arization induced by massive quarks has recently been 
computed in four-loop approximatioiP^ more precisely: its first derivative at = has been 
evaluated, which corresponds to the lowest moment of the familiar i?-ratio. With the help of the 
traditional integration-by-parts method in combination with Laporta's algorithm^ ^ a ll fo ur-loop 
integrals were reduced to a small set of master integrals which were taken from Refs.^^^. Based 
on these developments a new determination of the quark masses has been performed in Ref.^. 

The extraction of mg from low moments of the cross section a{e^e~ — > QQ) exploits its 
sharp rise close to the threshold for open charm and bottom production and the importance of 
the contributions from the narrow quarkonium resonances. By evaluating the moments 
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with low values of n, the long distance contributions are averaged out and A4n involves short 
distance physics only, with a characteristic scale of order -E'threshoid = 2mQ. Through dispersion 
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Figure 1: R{s) in the charm threshold region. 



relations the moments are directly related to derivatives of the vacuum polarization function at 

(?2 = 0, 
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which can be evaluated in perturbative QCD (pQCD). 

The narrow charmonium resonances J/^, ^'(25') and the higher excitations will obviously 
contribute to the moments. Open charm production exhibits a sharp rise, nearly like a step 
function. Beyond the ^'(3770)-resonance a few oscillations are observed which quickly level 
out into a fairly flat energy dependence. Around and above approximately 5 GeV the cross 
section is well ap proximated by pQCD and, furthermore, mass terms can be considered as small 
correction s ^ ' ^^ [ The sensitivity to uiq is, therefore, concentrated on the small region from J/^ 
up to approximately 5 GeV. 

We therefore distinguish three energy regions: First, the region of the narrow resonances 
J/^ and ^'(25'), second, the "charm threshold region" starting from the D-meson threshold at 
3.73 GeV up to approximately 5 GeV, where the cross section exhibits rapid variations and, third, 
the continuum region where pQCD and local duality are expected to giv e reliable predictions. For 
the threshold region we use the data from the BES collaboration^l^^I^, shown in Fig. [1] together 
with data from MD-l^I^l a,nd CLEO^^. Evidently pQCD provides an excellent description of 
all the data in the continuum region. The description of the perturbative continuum includes 
the complete mass dependence up to 0{al) plus the dominant mass dependent 0{ag) terms^ 
which were used to extrapolate Ruds from the region below charm threshold up to 4.8 GeV. 

In its domain of analyticity Ilc{q'^) can be cast into the form 
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where rric = 'mc{fj,) is the MS charm quark mass at the scale fi. The perturbative series for the 
coefficients Cn in order was evaluated in Ref.^^,the four-loop contributions to Cq and Ci in 



Refs.'^^. The coefficients depend on and on the charm quark mass through logarithms of 
the form = ln(m^(ju)/^^). Combining Eqs. ([T]), ([2]) and ([3]), the charm quark mass can be 
obtained: 
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Figure 2: Comparison of recent determinations of mi,(mi,) 



In the charm threshold region (which includes \I'(3770)) we have to identify the contribution 
from the charm quark, i.e. we have to subtract the parts arising from the light u, d and s quark. 
In the continuum region above \/s = 4.8 GeV data are sparse and imprecise. On the other 
hand, pQCD provides reliable predictions for R{s). Thus in this region we replace data by the 
theoretical prediction. 

We use the results for the moments to obtain in a first step mc(3 GeV). The moment with 
n = 1 is least sensitive to non-perturbative contributions from condensates, to the Coulombic 
higher order effects, the variation of ^ and the parametric dependence. We therefore adopt 

mc(3 GeV) = 0.986(13) GeV , (5) 

as our final result. Transforming this to the s cale- in variant mass mcirric) one finds mdmc) = 
1.286(13) GeV. Using the three-loop relation ^l^^lSl between pole- and MS-mass this leads to 

^^{3-loop) ^ ^ ggg 

The same approach is also applicable to the determination of m^. Just as in the charm case, 
a remarkable consistency and stability is observed. For n = 1 the error is dominated by the 
experimental input. For n = 3 we obtain ±0.010 from the experimental input, ±0.014 from 
and ±0.006 from the variation of /x. The three results based on n = 1, 2 and 3 are of comparable 
precision. The relative size of the contributions from the threshold and the continuum region 
decreases for the moments n = 2 and 3. On the other hand, the theory uncertainty estimated 
from the variation of fi and the unknown four-loop contribution is still acceptable. Therefore 
the result from n = 2 is taken as the final answer, 

mb(10 GeV) = 3.609(25) GeV, (6) 

corresponding to mbirrib) = 4.164(25) GeV and a pole mass of Adj^^ = 4.800 GeV. A 

comparison with a few selected determinations is shown in Fig. [5J 



For various applications, either related to Z-boson decays or in connections to Grand Unified 
Theories (GUTs), the values of at Mz and mt{mt) = 161.8 ± 2.0 GeV (as derived from 



rribiMz) = 2.834 ± 0.019 ± 0.017 GeV , mfe(161.8) = 2.703 ± 0.018 ± 0.019 GeV . (7) 

The first error reflects the combined error from Eq. ^ and the second one the uncertainty 
due to as- The ratio mt{mt) / mb{mt) = 59.8 ± 1.3 should be a useful input for Grand Unified 
Theories. 
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